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This is a reply report tc lQ- We studied heavy element production in the high baryon density 
region in the early universe[3|7 However it is claimed by [l[ that a small scale but high baryon density 
region contradicts observations for the light element abundance or, in order not to contradict to the 
observations the high density region must be so small that it cannot affect the present heavy element 
abundance. 

In this paper we study big bang nucleosynthesis in the high baryon density region and show that 
in certain parameter spaces it is possible to produce enough amount of the heavy element without 
O |. contradiction to cosmic microwave background and light element observations. 



PACS numbers: 26.35.+C, 98.80.Ft, 13.60.Rj 



I. INTRODUCTION 



In a standard scenario, big bang nucleosynthesis (BBN) can produce only light elements, up to 7 Li, and all heavy 
elements have been synthesized in stars. However, many phase transitions in the early universe could have printed 
(f") | their trace in a non-standard way. For example, some baryogenesis models0 predict very high baryon density islands 
\& • in ordinary low density backgrounds. 

In the previous paper we studied heavy element production in inhomogeneous BBN from this point of view. 1 
t^J- . However we limited ourselves to the heavy element abundance and did not discuss about the light element abundance 
and consistency with observations. This is because we assumed that high baryon density region is very local and 
r*"** , do not affect the global light element abundance. In Rauscher pointed out that in order not to contradict 
to observations, the high baryon density region must be very small and cannot affect the present heavy element 
abundance. In this paper, we show that there is a parameter region in which the heavy element can be produced 
enough to affect observation while keeping the light element abundance consistent with observations. We consider 
that the disagreement between Rauscher's opinion and our opinion comes from two points. One is that we are looking 
at some parameter regions in which neutrons in high baryon density do not diffuse so much as to cause disaster in 
standard BBN. We would like to emphasize this point. The other is that the relevant quantity is not the spatial size 
of the high baryon density region but the amount of baryon in high density regions. 

We will discuss the following issues: In section [H] we discuss the light element abundance in the homogeneous high 
baryon density region and after mixing the high and the low baryon density region. In section ITTTl we study the heavy 
element(Ru,Mo) abundance in high and averaged baryon density and show that heavy elements can be produced 
enough without contradicting the light element observation. In section [TVl we briefly comment on the diffusion scale 
of the high baryon density region. 



1 For previous works on the inhomogeneous big bang nucleosynthesis, see 4, 5]. Heavy elements production is also mentioned in |(J. 
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II. LIGHT ELEMENT ABUNDANCE 



A. Homogeneous BBN 

We calculate homogeneous BBN with various values of ^(baryon photon ratio). In Table HI and ITTl we show the 
numerical result of the mass fraction and the number fraction of each light clement for 77 = 10" 3 and 3.162 x 10 -10 . 





T) = 10" 


•i 




name 


mass fraction 


number fraction 


H 


5.814 x 10" 1 


8.475 x 10 


-1 


4 He 


4.185 x 10" 1 


1.525 x 10" 


-1 


3 He 


4.842 x 1(T 13 


1.614 x 10" 


-13 


7 Li+ 7 Be 


1.559 x 10" 12 


2.227 x 10" 


-13 


D 


1.577 x 1CT 22 


7.883 x 10" 
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Table. I: The mass and the number fractions of light elements for the homogeneous BBN with 77 = 10 



77 = 3.162 x 10" 



name 


mass fraction 


number fraction 


H 


7.58 x 10" 1 


9.26 x 10" 1 


4 He 


2.419 x 10" 1 


7.39 x 10" 2 


3 He 


4.299 x 10" 5 


1.433 x 10" 5 


7 Li + 7 Be 


8.239 x 10" 10 


1.177 x 10" 10 


D 


1.345 x 10" 4 


6.723 x 10" 5 



Table. II: The mass and the number fractions of light elements for the homogeneous BBN with 77 = 3.162 x 10 



As baryon density becomes higher, more protons and neutrons are bounded to form 4 He. At r\ = 10~ 3 , most of the 
final product of 7 Li comes from 7 Be which decays to 7 Li after BBN. Details on light element production for various 77 
can also be found in Q- In this paper we almost concentrate on a case in which the high baryon density region has 
X] = 10" 3 . We expect that compared to 77 > 10" 3 , the profile of the abundance for 77 = 10" 3 is more different from 
standard BBN because most of the light element abundances change monotonically with respect to 77 and if this case 
does not contradict to observations, other cases would also be consistent. Briefly, the amount of H decreases and 4 He 
increases monotonically as 77 become larger. The number fraction of D is less than 10~ 20 for 77 greater than 10" 7 . For 
3 He, the number fraction drastically decreases around 7/ = 10" 4 down to O(10 -13 ), and for 7 Li, the number fraction 
increases until 77 = 10" 6 and drastically decreases for a larger value of 77. In the following sections, we will see that 
this non-standard setup does not strongly contradict to the observations. For simplicity we ignore the diffusion effect 
before and during BBN, and after BBN both high and low baryon density regions are completely mixed. Detailed 
analysis such as the case in which the high baryon density region doesn't completely mixed, or taking into account 
diffusion effects are left for future work. 



B. Parameters and Basic equations 



In this section, we summarize the relations among parameters. 

Notations : n, n , n L are averaged, high, and low baryon number density. / , f L are the volume fractions of the 
high and the low baryon density region, yi, yf, are the mass fractions of each element (i) in averaged-, high- and 
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low-density regions. The basic relations are 

1 (1) 

(2) 

Din. (3) 
Under the assumption that the temperature of the universe is homogeneous, the above equation can be written as 

f H V H + / V = V (4) 

yff H v H + Vif L v L = ViV (5) 

where r\ = -^-,r] H ' L — " n Conventional parameters for inhomogeneous BBN are 77, / and density ratio R = '^ r . Here 

we use a different combination of parameters. Relevant values for the abundance analysis are products f H L x rj H ' L 
and r] H,L . fy' L x r] H ' L determines the amount of baryon from high- and low- density regions. rj H ' L determines the 
mass fraction of each species of nuclei. For convenience, we write the ratio of baryon number contribution from high 
density region as a, i.e., j H rj H : f L i] L — a : (1 — a). There are 5 parameters(n ff:i , n and f H ' L ) and 2 constraints 
(Eq.([T]) and Eq.([2])). We calculate the light element abundance for various values of rj H ' L . r\ can also take any value, 
but in order not to contradict observational constraints, we choose 77 from 3.162 x 10" 10 to 1CT 9 . a is determined 
by Eq(j4|. The aim of the analysis in this section is not to find parameter regions which precisely agree with the 
observational light element abundance and 77 from CMB. Our model is too simple to determine the constraints to 
parameters. For example, we completely ignore the diffusion effect before and during BBN. Instead we see that at 
least our analysis in previous paper is physically reasonable. 



f H + f L = 
f H n« + f L n L = 

y"f H n H + ytf L n L = 



C. Theoretical predictions and observations of light elements 

We consider the cases of r\ H = 10~ 3 and i] L = 3.162 x 10 -10 . The mass fractions of and H and 3 He in the 
high density region are 0.5814 and 4.842 x 10~ 13 , respectively, while those in the low density region are 0.758 and 
4.299 x 10~ 5 . From Eq.©, we have 

f H v H y"u e + f L v L ys Ue = wne (6) 

4.842 x 10" 13 x a + 4.299 x 10~ 5 x (1 - a) = y 3Rc (7) 

f H v H yg + f L v L vk = WH (8) 

0.5814 x a + 0.758 x (1 - a) = y H . (9) 
We can calculate an averaged value of the abundance ratio of 3 Hc to H as 

3 He _ 1 4.842 x 10~ 13 x a + 4.299 x 10~ 5 x (1 - a) 



where a is related to rj as 



H 3 0.5814 x a + 0. 758 x (1 - a) 



i] H rj — r/ L 
77 rj H — r\ L 

10- 3 77 - 3.162 x 10~ 10 

77 10- 3 - 3.162 x 10" 10 
7; - 3.162 x 10" 10 



(11) 
(12) 

(13) 



Here a varies from to 0.9 for reasonable values of 77, or 3.162 x 10 10 — 10 9 . Similarly, for 77^ = 10 3 the number 
fractions are 

,D, 1 1.577 x 10~ 22 x a + 1.345 x 10~ 4 x (1 - a) 

(—) = (14) 

V H ; 2 0.5814 x a + 0.758 x (1 -a) v ; 



1 .Oe-04 
8.0e-05 
6.0e-05 - 
4.0e-05 
2.0e-05 - 
O.Oe+00 - 



5e-10 7e-10 1e 



Fig. 1: Averaged ratio of D to H,(D/H) vs 77 
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Fig. 2: Same as Fig[T]but for ( 3 He/H) 
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Fig. 3: Same as Fig[T]but for ( 7 Li/H) 



5 



7 Li _ 1 1.559 x 1CT 12 x a + 8.239 x 1CT 10 x (1 - a) 

(~H~> ~ 7 0.5814 x a + 0.758 x (1 - a) ' ^ ' 

Fig lll2l and [51 represent the averaged abundance ratio, (D/H), ( 3 He/H) and ( 7 Li/H) respectively. 

We can see that the light element abundance is the same order around r\ ~ 5 x 10~ 10 — 10~ 9 as observations 

ME [II EH Bill. 

& obs = (1.5-6.7) x 10- 5 (16) 
ri 

(-^Us = (0.59-4.1) x 10- 10 . (17) 

We do not discuss detail about diffusion here. But at least above result suggest that our analysis is not beside the 
point. 



III. THEORETICAL PREDICTIONS AND OBSERVATIONS OF HEAVY ELEMENTS ( 92 94 MO, 96 98 RU) 

The same analysis can be applied for heavy elements such as 92 Mo, 94 Mo, 96 Ru and 98 Ru. We are interested in 
these elements because in many models of supernovae nucleosynthesis, these p-nuclei are less produced. We will see 
that some amount of these heavy elements can be synthesized in BBN. 



Tj = 10" J 

name mass fraction 



H 


5.814 x 10" 


i 


4 He 


4.185 x 10" 


i 


92 Mo 


1.835 x 10" 


5 


94 Mo 


4.1145 x 10 


-6 


96 Ru 


1.0789 x 10 


-5 


98 Ru 


1.0362 x 10 


-5 



Table. Ill: The mass fractions of nuclei for homogeneous BBN with r\ — 10 



From Table \IU\ we can derive the expected value of these elements. 



92 Mo, 1 1.835 x 10" 5 x a 

( tt ) - 777 n eoi a Z 1 n 7 C o ~ Ti \ ( 18 ) 



H ' 92 0.5814 x a + 0.758 x (1 - a) 

4 Mo _ 1 4.1145 x 10" 6 x a 
~H~' ~ 94 0.5814 x a + 0.758 x (1 - a) 

6 Ru _ 1 1.0789 x 10" 5 x a 

~H~' ~ 96 0.5814 x a + 0.758 x (1 - a) 

^Ru _ 1 1.0362 x 10~ 5 x a 

B~' ~ 98 0.5814 x a + 0.758 x (1 - a) ' 



94 Mo N 1 4.1145 x 10~ 6 x a , x 

(^) = 771 7777777777 I 77777T7777 ^ (19) 



^ 96 Ru ^ _ 1 1.0789 x IP" 5 x a ^ 



98 Ru 1 1.0362 x 10" 5 x a 

> xj ' no n coi A ^, „ i n 7t;o \, ft „\ ' \ J 



We plot expected value of these quantities in FigHJ 

These values should be compared with the solar abundance ( Table II VP 17]. 
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Fig. 4: ( 92 Mo/H),( 94 Mo/H),( 96 Ru/H) and ( 98 Ru/H) vs 77. Red, green, blue and pink lines represent the ratio 
( 92 Mo/H),( 94 Mo/H),( 96 Ru/H),( 98 Ru/H) respectively. 



name 


number fraction 


ratio to H 




H 


7.057280 x 10 


-1 


1 




92 Mo 


8.796560 x 10" 


-10 


1.2465 x 10" 


-9 


94 Mo 


5.611420 x 10" 


-10 


7.9512 x 10" 


10 


96 Ru 


2.501160 x 10" 


-10 


3.5441 x 10" 


10 


98 Ru 


8.676150 x 10" 


-11 


1.2294 x 10" 


10 



Table. IV: The abundances of 92 ' 94 Mo and 96 ' 98 Ru in the solar system[T3] 




Compared those observational values with FigHJ it is clear that the heavy element produced in BBN can affect the 
solar abundance heavy element. Some of them are produced too much. But this is not a problem of the previous work 
[H, because we assumed that high density regions are very small and do not disturb standard BBN. The analysis here 
suggest that even if we assume the density fluctuations are completely mixed, heavy element can have enough affect 
to the solar abundance. 



IV. DIFFUSION DURING BBN 



In the previous analysis, we assumed that the diffusion effect can be ignored during BBN and both high density 
regions and low density regions are completely mixed after BBN. In this section, we determine the scale of high baryon 
density island in which the diffusion effect during BBN is very small enough and our assumption is valid. We do not 
discuss the diffusion after BBN here. 

A detail analysis of the comoving diffusion distance of the baryon, the neutron and the proton is in (l6| . From Fig.l 
in [l6| . in order to safely ignore the diffusion effect, it is necessary for the high baryon density island to be much larger 
than 10 5 cm at T=0.1MeV(l.l x 10 9 K). Notice that T oc where A is a scale factor. For scale d now corresponds to 
(f/(4.0 x 10 8 ) at BBN epoch. Present galaxy scale is O(10 20 )cm, which corresponds to O(10 12 )cm >> 10 5 cm at BBN 
epoch. 



temperature and scale 


temperature 


scale 


1.1 x 10 9 K (BBN) 


d 


3000K (decouple) 


3.7 x 10 6 x d 


2.725K (now) 


4.0 x 10 8 x d 



Table. V: Relation between temperature and scale 
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The maximum angular resolution of CMB is l ma x ^2000. The size of universe is ~ 5000Mpc. In order not to 
contradict to CMB observation, the fluctuation of baryon density must be less than ~ 16Mpc now. This corresponds 
to 10 17 cm at BBN. 

Since the density fluctuation size in Dolgov and Silk's model[2j is a free parameter, the above brief estimation sug- 
gests that we can take the island size large enough to ignore the diffusion effect without contradicting to observations, 
i.e., the reasonable size of 10 5 cm — 10 17 cm at the BBN epoch. We can choose distances between high density islands 
so that we obtain a suitable value of /. 

V. SUMMARY 

In this paper, we studied the relation between the heavy element production in high baryon density regions during 
BBN and the light element observation. By averaging the light element abundances in the high and the low density 
regions we showed that it is possible to produce a relevant amount of heavy element without contradicting to observa- 
tions. However we should stress that in this paper we restricted ourselves to some parameter regions where neutrons 
in high baryon density regions do not destroy the standard BBN. So our setup is different from the conventional 
inhomogeneous BBN studies. We also studied the size of the density fluctuation to show that there is a parameter 
region in which the neutron diffusion is negligible and which is much smaller than CMB observation scale. It is 
worthwhile to investigate further how the produced heavy elements can be related to the detailed observations. 
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